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Abstract 

The 5't/(3)-violating decays ^^^X, where X = or rj and ^ = J/tp or T have been 

recently proposed as a means of probing the hght quark masses beyond leading order in 
chiral perturbation theory. We argue that this analysis is incorrect, even in the heavy quark 
limit. We show that these decays are governed by an infinite number of matrix elements 
which are not suppressed by any small parameter, and which cannot be computed with our 
present understanding of QCD. Furthermore, for sufficiently heavy quarks, we show that 
the decay amplitudes can be organized into a twist expansion, and that the contributions 
considered in the above proposal are subleading in this expansion. We also explain how 
these decays nonetheless give a constraint on the light quark masses valid at leading order 
in the chiral expansion. The decays ^r)'y and ^"^^ ^ ^^^tttv also have contributions 
from infinitely many operators, contrary to claims in the literature. 
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1. Introduction 



Recently, there has been a renewed interest in the flavor SU (3)-violating decays 

where ^ = J/ijj or T, as a probe for the light quark masses beyond leading order in chiral 
perturbation theory [1] [2] . This work is based on ref. [3] , where it is argued that the QCD 
multipole expansion [4] [5] can be used to show that 

^($25^^isj^) ^ (X\tTG'"'G^,,\0), (2) 

where X = tt'' or r]. The basic idea in ref. [3] is that the $ is spatially small compared 
to the wavelength of gluons which arc important in the decay process, and that the decay 
can be computed in terms of the matrix elements of local gluonic operators. This result is 
supposed to be valid for sufficiently heavy quark masses. The hadronic matrix element in 
eq. (2) can then be computed using chiral perturbation theory [6] [2]. 

In this paper, we show that eq. (2) is not the leading contribution to the decays in 
eq. (1): there are an infinite number of operators which contribute whose matrix elements 
are not suppressed by any small parameter. This result is a simple consequence of the 
non-locality of the transition amplitude in time. Furthermore, in the heavy quark limit, 
the decay amplitude can be organized into a twist expansion, and the contribution leading 
to eq. (2) can be shown to be subleading in this expansion. This invalidates the claim 
of ref. [2] that the decays in eq. (1) can be used to constrain the light quark masses at 
next-to-leading order in the chiral expansion. Despite this negative result, we show that 
the decays (1) can be used to derive a constraint on the light quark masses valid at leading 
order in the chiral expansion, but not at higher orders. This distinction is important, since 
the corrections to lowest-order results are substantial, and may even allow m-u = 0, giving 
an economical solution of the strong CP problem [7]. 

A similar analysis can be used to show that the decays 

$2S^^is^^ and $^^^7/7 (3) 

also get contributions from inflnitely many operators, contrary to the claims of refs. [3] 
and [6], respectively. 

2. Quantum Mechanics Analysis 

In this section, we adopt the assumptions and methodology used in ref. [3] to describe 
the decay in eq. (1). We show that these assumptions lead to the conclusion that these 
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decays are governed by an infinite number of matrix elements which are unsuppressed by 
any small parameter. In the next section, we will argue that this conclusion is much more 
general than the argument of this section, and that the contributions considered in this 
section are in fact not the dominant ones. 

We will treat the $ as a non-relativistic quantum-mechanical bound state. We begin 
by making the standard assumption that the decay in eq. (1) can be viewed as a two- 
step process: First the decays via the emission of two gluons, and then the gluons 
"hadronize" to form the light pseudoscalar in the final state. The wavelength of the gluons 
involved is large compared to the spatial size of the and so gluon emission can be 
described in terms of a local perturbation hamiltonian involving the gluon fields. 

In ref. [3], it is assumed that the leading contribution to the decays in eq. (1) comes 
from interference between the electric dipole interaction 

SHi = -gf-E, (4) 



and the the spin-dependent interaction 

5H, = --^{r-D)iS-B), (5) 

which arises at higher order in the nonrelativistic expansion. Here, r is the separation 
between the quark and the antiquark, S is the total spin operator, and 

E = -VA^^Ta, B = Vx AaTa, (6) 



where 

T„=i(t„®i + i®t:). (7) 

The matrix element is then given in second-order perturbation theory by 

A oc (/I'^^i— -TTr^TTTTTnT ^^^'^^ + 2), (8) 

where is the energy of the state is the hamiltonian of the QQ pair in the 

QCD-generated potential and Hg is the hamiltonian for free gluons. Inserting a complete 
set of energy eigenstates, we obtain 

A oc ^(/l^ifi _ „ _ ^ ——\a){a\6H2\i) + (1 - 2), (9) 
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where uia. is the energy carried away by the emission of the first gluon. We can expand 
this to obtain 



A^YJZWa^ -VzO+)n+i Wm^a5H,\i) + (1 - 2) 

n a ^ * QQ ' 



(10) 



where the time derivative is understood to act on the gluon fields in ^H^. We therefore 
obtain 

^ ^ 2mQ^^ ' ^(e,-iyQQ+zO+)-+i (11) 

X {X\Ej{iDorDkBe + DjBe{iDorEk\0), 

where we have included additional terms mandated by gauge invariance by substituting 
the gauge-covariant derivative Dq for the time derivative Oq. Similar results are derived 
using field-theoretic arguments in ref. [5]. 

If we expand this result in powers of Hqqi eq. (11) has the form of an operator 
product expansion, with the suppression scale of the higher-dimension operators given 
by the spacing between excited energy levels A. If we keep only the first term in this 
expansion, we recover the result eq. (2) of ref. [3]. However, the kinematic energy scale for 
these decays is M2S — Mis, which is of order A, so there is no small parameter suppressing 
the higher order terms in the series in eq. (11), even in the heavy-quark limit. This is the 
central point of this paper. 

To show explicitly that there is no suppression of higher-order terms in eq. (11), we 
follow ref. [3] and assume that the spatial and spin parts of the $ wavefunction factorize. 
We can then write eq. (11) as 

^ ^ 2mQ^A-+3^ (12) 

X {X\Eu{iDQYDkBj + DuB^{iDoTEk\Q), 

where is a quantum-mechanical matrix element, and and are polarization 
vectors for the initial and final quarkonium states. In order to estimate the size of the 
hadronic matrix element, wc write the operator in relativistic notation using the four- 
velocity V of the decaying particle: 
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where 



Here, we use the notation (■ • •) to denote the symmetrization of indices. The hadronic 
matrix elements can be written 

(^l^^i;U...AjO) = [g^,^p^, . . -vx^^v, - - P)] F^-^ ^^^^ 

+ \g{fxv9\rMV\z ■ ■ •PA„)Pp -{v^ p)] + • • • , 

where p is the momentum of X. We have expHcitly factored out powers of the hadronic 
scale A so that all the form factors Fj have the same dimension. Since the form factors 
depend only on jp" = m^, they are constants determined by low-energy QCD dynamics, 
and are all expected to be of the same order of magnitude. (Note that since the decay 
violates SU{3), the form factors Fj must be proportional to quark masses.) 

Substituting this into eq. (13), we obtain 

^($25 ^ $i^X) oc ^ A^ v.e^'^'^^l'^fp, {p ■ vT Ft^ + ■■■. (16) 

Since p-v ^ A, we see that the contributions from the higher order terms in this expansion 
are unsuppressed, even in the heavy quark limit. 

Similar arguments which led to eq. (2) were used in ref. [3] to argue that 

^($2S ^ ^iStttt) oc (ttttI trG'^'^G^^lO), (17) 
and in ref. [6] to argue that 

A{^^s^rj^) oc (r7|trG^-G^,|0). (18) 

The considerations above can be used to show that these decays also have unsuppressed 
contributions from an infinite number of operators. 



3. Effective Field Theory Analysis 

The conclusions of the previous section are not at all unexpected if we consider these 
decays in an effective field theory language. If we consider interactions of the state 
with gluons with energies <^ A, then we can write an effective field theory in which the 

is treated as a heavy particle. Higher-dimension terms in this effective lagrangian are 
suppressed by powers of A, the energy required to excite the next excited state. In order to 
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consider decays such as those in eq. (1), we must write an effective theory containing fields 
for both and states interacting with gluons. This is not a very useful effective 
theory, since the higher-dimension terms in this lagrangian are again suppressed by powers 
of A, which is also the kinematic scale for the decays we are interested in. Therefore, we 
expect that an infinite number of gluonic operators will contribute at the same level to the 
decay. This situation arises because we are trying to describe a process which is nonlocal 
(in time) on a scale ~ 1/ A by local operators. 

In an effective field theory for $ radial transitions, we could immediately write down 
contributions to the decay such as eq. (13) and reach the conclusions of the previous 
section. But we can also write terms such as 

= ^v^.^^'^'K'^l' tr(Gf,l?A. • • • D^^D,GfS)a)v''v^v>^^ • • • + • • • . (19) 

There is no symmetry which forbids these terms, so we expect them to be generated at 
some level. In fact, at higher orders in the non-relativistic expansion, factorization of the 
spin and spatial parts of the ^ wavefunction assumed in eq. (12) breaks down, generating 
the contributions in eq. (19). These terms give rise to contributions to the decay amplitude 
proportional to 

{X\tr{Gl^Dx, ■ ■■Dx^DpGf3),)\0) = PaPpPx, ■ ■ -Px^PpGi + ■■■, (20) 

where the form factor has been normalized to have the same dimension as those in eq. 
(15). In the heavy-quark limit, the contribution of this term to the decay amplitude is 
enhanced relative to eq. (16) by a factor of 

up to logarithms of mg, where A ~ 1 GeV is a hadronic scale. 
In general, for terms in the effective theory of the form 

^eff = E 1^ v^e^^'^K'^fO^;^ , (22) 

n 

where O^"^ is a QCD operator, it is easy to see that the leading contributions will come 
from operators with lowest twist (= dimension— spin). The operator in eq. (15) has twist 4, 
while the operator in eq. (20) has twist 2. Thus, the contributions to the decay amplitude 
in eq. (12) are not among the leading ones in the heavy-quark limit. 

The relevance of the twist expansion to the real world is unclear, since the c and b 
quarks are not sufficiently heavy to trust heavy-quark results without reservation, and in 
any case the hadronic matrix elements of the operators in question are unknown. 
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4. Determination of the Light Quark Masses 

Our main motivation for this work is the fact that ref. [2] argued that the decays in 
eq. (1) can be used to give a constraint on the hght quark masses beyond leading order in 
the chiral expansion. We have argued that these decays are not dominated by the matrix 
elements in eq. (2), so it may seem mysterious that the constraint obtained in ref. [2] is 
consistent with other determinations of the light quark mass ratios. However, it is easy to 
see that the decays (1) give a constraint on the light quark masses valid at leading order 
in the chiral expansion. 

The point is that the decay is forbidden in the limit of exact vector SU{3). However, 
SU (3) is broken explicitly by the quark masses, so the leading contribution to the amplitude 
can be written 

^($2S ^ oc tY{MqTx), (23) 

where Mq is the light quark mass matrix and Tx is the SU{3) generator corresponding to 
the pseudoscalar in the final state. (We can use chiral invariance to show that a similar 
term does not exist due to electromagnetic breaking of SU{3).) This form of the amplitude 
immediately gives the leading order constraint on the light quark masses used in ref. [2]. 
Similar arguments to those given above can be found in ref. [1]. 

5. Conclusions 

We have shown that the decays in eq. (1) cannot be used to constrain the light 
quark masses beyond leading order. These decays are governed by an infinite number of 
operators which are unsuppressed by any small parameter even in the heavy quark limit. 
Furthermore, for sufficiently heavy quark masses, the leading contributions to the decay 
come from operators of twist 2, while the contribution discussed in ref. [3] has twist 4. 

The use of the decays in eq. (1) to determine the light quark masses has also been 
criticized in ref. [8]. However, this paper assumes that the leading contribution to the 
decay comes from the first term of eq. (12), and criticizes the quantitative accuracy of the 
heavy quark and 1/A^coior expansion used in ref. [3]. 
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